Hydrodynamic modes in a trapped Bose gas above the Bose-Einstein transition 
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We discuss the collective modes of a trapped Bose gas in the hydrodynamic regime where atomic 
collisions ensure local thermal equilibrium for the distribution function. Starting from the conser- 
vation laws, in the linearized limit we derive a closed equation for the velocity fluctuations in a 
trapped Bose gas above the Bose-Einstein transition temperature. Explicit solutions for a parabolic 
trap are given. We find that the surface modes have the same dispersion relation as the one recently 
obtained by Stringari for the oscillations of the condensate at T — within the Thomas-Fermi 
approximation. Results are also given for the monopole "breathing" mode as well as for the m — 
excitations which result from the coupling of the monopole and quadrupole modes in an anisotropic 
parabolic well. 



PACS numbers: 03.75.Fi, 05.30.Jp, 67.40.Db 

It is well known in discussions of uniform quantum flu- 
ids that there are two kinds of collective modes Q . There 
are collisionless (zero sound) modes arising from dynamic 
self-consistent mean fields. Recently there have been sev- 
eral theoretical studies [|],[| of the collective modes of a 
trapped Bose gas at T = 0, where all the atoms are in the 
condensate These are all based on the Gross-Pitaevskii 
description corresponding to the time-dependent Hartree 
theory of the Bose condensate, which can be also gen- 
eralized to finite temperature Recently the lowest 
frequency oscillations of a trapped Bose condensate have 
been measured and found to be in excellent agree- 
ment with the calculated T — Bogoliubov excitation 
frequencies [§,1). 

In addition, there can be hydrodynamic modes (such 
as first and second sound in Bose superfluids) which arise 
when collisions are sufficiently strong to ensure local ther- 
modynamic equilibrium. This letter is concerned with 
the hydrodynamic modes of a trapped Bose gas. We de- 
rive a closed equation for the velocity fluctuations in a 
gas which is not Bose-condensed (T > Tbec)- We ex- 
hibit exact solutions of this equation for a parabolic well 
corresponding to surface, monopole and m = coupled 
monopole-quadrupole modes. The frequencies of these 
normal modes are shown to be the same for a degener- 
ate Bose gas as for a classical gas. The surface modes 
have frequencies identical to those found by Stringari || 
for a T = Bose condensate, in the Thomas-Fermi ap- 
proximation ||, while the monopole mode has the same 
frequency as in a non-interacting gas. 

The hydrodynamic regime is achieved when the colli- 
sions! frequency u c — l/r c is much larger than the fre- 
quency bj of the collective mode {ojt c <C 1). For ui ~ too 
(the trap frequency), this condition is very similar to re- 
quiring that the size of the system R be much larger than 
the collisional mean free path I ~ 1/na, where a — 8ira 2 
is the s-wave cross section and n is the density Q. At 
T ~ 2Tbec, these conditions seem to be well satisfied 



in a recent MIT experiment ||, while in the JILA ex- 
periment |(|] the system is too dilute and is in the non- 
interacting regime above Tbec- 

Kadanoff and Baym have given a very clear treat- 
ment of collision-dominated hydrodynamic modes in a 
uniform dilute classical gas which is easily extended to a 
degenerate Bose gas in an external potential Uo(r). The 
collisions ensure that the perturbed distribution function 
produced by a slowly varying external field SU (r, t) is al- 
ways given by the local thermodynamic equilibrium Bose 
distribution function 



/(p,r,i) = (< 



/3(r,t)[(p-mv(r,t)) 2 /2m-M(r,i)] 



(1) 



The perturbed distribution is described by the space- and 
time-dependent temperature T(r,t)) = [fes/3(r, t)] _1 , 
chemicalpotential fi(r,t), and local velocity v(r,t). In- 
serting (p into the conservation laws for the local particle 
number n(r, t), particle current j(r, t) and energy density 
e(r, t), one finds equations of motion which determine the 
values of T(r, t), fi(r, t) and v(r, t). (This procedure is 
described on p. 56-58 of Ref. |1(J for a uniform classical 
gas). We shall only be interested in small perturbations 
around the equilibrium state, 



T(r, t) = T + ST{r, t) 
H(r,t) = /Lt (r) +S(j,(T,t) 
v(r,t) = 5v(M) • 



(2) 



We note that in the equilibrium state of the non-uniform 
trapped gas, the temperature Tq is uniform and the local 
velocity vo(r) is zero. 

Working to first order in 8v(r, t), the conservation laws 
for the local perturbed quantities 



n(r, t) 
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reduce to [nO 



mno (r) 



dn(r, t) 
dt 



d6v(r,t) 
dt 



-V-[no(r)Jv(r,t)] 



(4) 



-[VP(r,t)+n(r,t)Vl7o(r)] (5) 



fls(r, f) 
9t 



-e (r)<5v(r, i) - n (r)5v(r, t) • VtT (r). 

(6) 

The terms proportional to V[/o(r) in (g) and @ show the 
explicit role of the trapping potential. The equilibrium 
density no(r) and pressure Po(r) must be consistent with 
the vanishing of the r.h.s. of (g) when <5v(r, t) = 0, 



VPo(r)+n (r)VC/ (r) =0. 



(7) 



Using in (H|), the integrals over p are exactly the same 
as for a uniform Bose gas in thermal equilibrium |TT| and 
one may easily verify that 



j(r,i) = no(r)<Sv(r,t) 



1 



P(r,t) 



1 
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(8) 
(9) 
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Here z(r, i) = e^ r ' t '^ r ' t ^ is the local thermodynamic 
equilibrium fugacity, the thermal de Broglie wavelength 
A = (2nh 2 /mk B T(r,t))^ 2 : and g n (z) =Y?L 1 z l /l n are 
the well-known Bose-Einstein functions |llf| . 

While collisions are crucial to enforce the local thermo- 
dynamic equilibrium solution given by ([l]), for a dilute 
gas one can, to lowest order, ignore the interactions in 
evaluating the hydrodynamic equations of motion from 
the conservation laws, as we have done in the preced- 
ing analysis. The static local equilibrium values of the 
thermodynamic functions are given by by set- 



ting <Sv(r, t) = 0, T(r, t) = T and z = z 



= "o(r)/fc B T 0) 

where /Lto(r) = \i — Uo(r) and \i is the chemical po- 
tential. We note that the temperature To is constant 
throughout the trap. These results correspond to the 
well-known semi-classical approximation |12], which is 
valid when the discrete energy level spacing of the trap- 
ping potential is much less than UbTq. The resulting 
equilibrium values are consistent with the relation given 
in (^), as easily verified using the well-known identity 
dg n {z)/dz = g n _x(z)/z. 

One can use to derive a closed equation for 

the velocity fluctuation Jv(r, t). Using (|l0|), (||) can be 
rewritten as 



= -|V • [P (r)Mr,i)] - ln (v)Sv(r,t) ■ VU (r). 

(11) 

Taking the time derivative of fl|) and using ( pTT| ) and (^) 

gives 



mn a (r) 



d 2 Sv 
dt 2 



^V[V ■ (P (r)5v)] + ^V[n (r)6v ■ Vf/ (r)] 



V-[n (r)5v]VC/ (r) 



(12) 



where 5v = Sv(r,t). Using (R) to rewrite V£/o(r) and 
the fact that Vno(r) oc V£/o(r), it is straightforward to 
reduce (12) to 



5v _ 5P (r) 



dt 2 3 n (r) 



V[V • 5v] - V[<5v • Vt/ (r)] 



-[V-<5v]VC/ (r) 



(13) 



This is the key result of this letter. We note that the 
equilibrium properties of the non-uniform Bose gas only 
enter into ( |l3| ) through the ratio 



Po(r) 
n (r) 



55/2 (^o) 
03/2 (*<>) 



(14) 
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and 
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where the local fugacity zo(r) has been already defined. 

In the absence of a trapping potential we have a uni- 
form gas and ([l^) has the plane wave solution <5v(r, t) = 
(5v w (k)e l ( k - r -^ with <Sv w (k) - k. This is a longitudinal 
hydrodynamic sound wave with the dispersion relation 
uj 2 = c 2 k 2 , where the sound velocity is given by 

c 2 _ J^_Po _ 5 k B T g 5 /2(zo) 
3m n 3 m g 3 / 2 {zo) 

In the classical limit, where zq — e At / fcBT ° 
duces to the well known Laplace result [jll| 
At Tbec, we have zq = 1 since ^ — 
c 2 = 0.51c 2 ,. 

We note that ( |l3|) is a vector equation, giving three 
coupled equations for the three components of the ve- 
locity 8v(r,t). In the presence of the external trap, the 
normal mode solutions of this hydrodynamic equation 
are not automatically irrotational and may contain rota- 
tional components. This is a consequence of the inhomo- 
geneous nature of the system. 

Special and important solutions of the form <5v(r, t) = 
(5v w (r)e _Mt can be obtained which are simultaneously 
irrotational and divergence free. For these normal modes, 
(O) reduces to 



then 



mw 2 <5v w (r) = -V[<Jv w (r) • VC/ ] 



(16) 



Since the first term in (|13|) makes no contribution, these 
solutions are the same in a classical gas as in a highly 
degenerate Bose gas just above Tbec- These zero- 
divergence solutions can be shown to involve isothermal 
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oscillations as follows. Using and (|T^), one finds 
(6 = k B T) 

dP 1 /5„ \ 96 9/i 

flT = e^U p °- n °^J -m +n °m (17) 

^ = ^U"°- 7o H^ +7o w (18) 



where 70 = A 3q gi/2(^o)- Combining these with (|l] 

and (Q), one can eliminate d[ijdt and solve for 80/ dt to 
obtain 



^e(r,t) = -He v-5v(r,t)- 



(19) 



This result, valid for all T > Tbec, shows that the tem- 
perature is constant for divergence-free modes. 

For an isotropic trap Uo(r) = ^mui 2 r 2 , a normal mode 
solution of ( |l~6| ) is given by 

6v u (r) ~ V[/F, m ] (20) 

with the dispersion relation u> = Viuio holding for I > 1. 
Using (||) one verifies that the associated density fluctu- 
ation $n(r, i) = (5rj w (r)e~ lwt is given by 



dr 



(21) 



In the classical limit, the static density profile reduces to 




n (r) = n (r = 0)e~ 



2 /2e 



(22) 



and the density fluctuation is then given by 8ti uj (y) ~ 
' ^m(6'j 0)"-o( r )- The latter vanishes at the origin and 
at large r, being peaked at r = {iQ^/muj 2 ) 1 / 2 . Such 
zero-divergence modes are usually referred to as surface 
modes. It is worth noting that the dispersion relation of 
these hydrodynamic surface modes is identical to the one 
found by Stringari || in an interacting Bose-condensed 
gas at T — when the equations of motion for the two- 
component order parameter take the form of the equa- 
tions of hydrodynamics Q . The equation for the velocity 
fluctuations in Ref. || can be shown to take the simple 
form 



d 2 Sv 



-V[5v ■ Vt/ (r)] + V[0u - U (r))V ■ 5v], (23) 



which reduces to ( |iq ) for zero-divergence modes. This ex- 
act correspondence is a surprise, but only at first sight. In 
fact, these surface modes are entirely driven by the exter- 
nal force and are insensitive to the form of the equation- 
of-state, which is quite different in the two cases [ p"3| . 
The dispersion relation of the surface modes discussed 
above differs from the excitations of an harmonic oscil- 
lator model in the absence of interactions, which are de- 
scribed by lu = £luq. 



We next discuss the compression mode solutions of 
(|l3|). The solution for the "breathing" monopole oscil- 
lation is obtained by setting (we assume an isotropic 
harmonic potential) (5v w (r) ~ r. Equation ( |l3|) then 
yields lo = 2loq and the density fluctuation takes the form 
<5n w (r) ~ (3— )n-o(r) in the classical limit [jl3). Be- 

cause V • <5v w (r) is a constant, this monopole solution is 
unaffected by the quantum terms associated with Bose 
statistics which enter ( |l3| ) only through the ratio Pq/tiq 
given by (|l4|). The analogous solution for the interact- 
ing Bose-condensed gas S at T = (see (p3|)) has the 
dispersion relation u> = v5^o, but in the non- interacting 
model the frequency is given by lo = 2loq. The fact that 
the monopole frequency coincides in the hydrodynamic 
regime and in the non-interacting model is a special fea- 
ture of a parabolic well. As first shown by Boltzmann 
PH , for an isotropic harmonic trap, the monopole oscil- 
lation of frequency 2wo in a classical gas can be shown 
to be an exact undamped solution of the full Boltzmann 
equation 15 1(| . 

Current Bose gas experiments [^]|| usually involve an 
anisotropic magnetic trap with axial symmetry 



C/o(r) = -m(^ S 2 + ^z 2 ) 



s 2 = x 2 + y 2 . (24) 



One easily finds that there are exact solutions of (|l^), 
with frequencies and associated density fluctuations (for 
a classical gas) given by 

w? )Tn=±I = lwi ; 6n u (r) ~ s l e ±u <t>n Q {v) (25) 



(26) 

where n (r) is the anisotropic version of (|||) using (|2J 
Identical mode frequencies were obtained by Stringari 
for an anisotropic condensate at T = 0, as described by 

In a deformed trap, the monopole mode is coupled to 
the quadrupole surface mode (£ = 2, m = 0). The decou- 
pling is easily carried out starting from (U3h and looking 
for irrotational solutions of the form (5vJ"[r) = Vx, with 
X = as 2 + (3z 2 . This gives two coupled equations for a 
and f3, which yield the dispersion relation 



Alj 2 ± ( 25^1 



16^-32^ 2 ) 1/2 ] 



(27) 



For an isotropic trap, this expression reduces to to = 
2wo and V2luo, corresponding to the monopole and 
quadrupole excitations, respectively. In a cigar-type con- 
figuration with lo z <C lo±, the lowest solution of d27]) has 
the frequency to 2 = ^^ 2 , very close to the analogous re- 
sult oj 2 — %uj 2 z for the oscillation of the Bose condensate 
at T = ||. This result is consistent with the recent 
data of the MIT group ||. They found essentially the 
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same normal mode frequencies ~ 30Hz at very low tem- 
peratures (T <C Tbec) and at T — 2 Tbec, the mode 
in the latter case being in the hydrodynamic region since 
v c ~ 4 x 10 3 Hz. In spite of the very similar frequen- 
cies, the damping of these collisionless and hydrodynamic 
modes is expected to be quite different Q|. 

In this paper, we have shown that the complete spec- 
trum of hydrodynamic normal modes in a trapped Bose 
gas above the BEC transition can be obtained from the 
solutions of the vector equation in ([l3]). We have used 
this equation to discuss the normal modes for which 
V • 8v uj (y) is constant, in which case the first term of 
the r.h.s. of ( |i~3| ) makes no contribution. This is the rea- 
son why the solutions are equally valid for a Bose gas just 
above Tbec as for a classical gas. We exhibited inter- 
esting similarities to the excitations of a Bose-condensed 
gas at T = and of a non-interacting gas. We hope that 
our results stimulate further experimental studies of the 
hydrodynamic modes in trapped Bose gases above the 
BEC transition. There are additional normal mode so- 
lutions of ( [Hf ) in which the first term involving Pq/uq 
does contribute. In this regard, we note that a plot of 
95/2(20)/ 53/2(^0) in dbj ) as a function of the distance 
r from the center of an isotropic parabolic well shows 
that this ratio is remarkably close to the classical value 
of unity down to about ITbec- At lower temperatures, 
it is only slightly less than unity at small values of r. 
Thus the solutions of ([l3]) which depend on the first term 
can be expected to be well approximated by the classical 
limit. These will be discussed elsewhere. 

In a separate publication, we will report on an exten- 
sion of the present work to the Bose-condensed region 
just below Tbec- This generalized version of ([|)-@ 
leads to the analogue of the two-fluid superfluid equa- 
tions jl7 18 1 for a trapped Bose gas, with a new degree of 
freedom associated with the superfluid component. The 
hydrodynamic modes (first and second sound) of a uni- 
form weakly interacting Bose gas have been discussed in 
Ref. U|. 

After completing this work, we received a preprint from 
Yu. Kagan, E.L. Surkovand G.V. Shlyapnikov. Their re- 
sults, based on the formalism of scaling transformations, 
agree with our expression (p7| ) for the quadrupole and 
monopole excitations. 

Stimulating discussions with Lev Pitaevskii and Eu- 
gene Zaremba are acknowledged. A.G. was supported by 
a grant from NSERC of Canada. 
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